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In modern societies, cultural change seems ceaseless. The flux of fashion is especially 
obvious for popular music. While much has been written about the origin and evolution 
of pop, most claims about its history are anecdotal rather than scientific in nature. To 
rectify this we investigate the US Billboard Hot 100 between 1960 and 2010. Using Music 
Information Retrieval (MIR) and text-mining tools we analyse the musical properties of 
^17,000 recordings that appeared in the charts and demonstrate quantitative trends 
in their harmonic and timbral properties. We then use these properties to produce an 
audio-based classification of musical styles and study the evolution of musical diversity 
and disparity, testing, and rejecting, several classical theories of cultural change. Finally, 
we investigate whether pop musical evolution has been gradual or punctuated. We show 
that, although pop music has evolved continuously, it did so with particular rapidity during 
three stylistic “revolutions” around 1964, 1983 and 1991. We conclude by discussing 
how our study points the way to a quantitative science of cultural change. 

Introduction 

The history of popular music has long been debated by philosophers, sociologists, journal¬ 
ists and pop stars [1-6]. Their accounts, though rich in vivid musical lore and aesthetic 
judgements, lack what scientists want; rigorous tests of clear hypotheses based on quanti¬ 
tative data and statistics. Economics-minded social scientists studying the history of music 
have done better, but they are less interested in music than the means by which it is mar¬ 
keted [7-14] . The contrast with evolutionary biology—a historical science rich in quantitative 
data and models—is striking; the more so since cultural and organismic variety are both con¬ 
sidered to be the result of modification-by-descent processes [15-18]. Indeed, linguists and 
archaeologists, studying the evolution of languages and material culture, commonly apply 
the same tools that evolutionary biologists do when studying the evolution of species [19-24]. 
Inspired by their example, here we investigate the “fossil record” of American popular music. 
We adopt a diachronic, historical approach to ask several general questions; Has the variety 
of popular music increased or decreased over time? Is evolutionary change in popular music 
continuous or discontinuous? If discontinuous, when did the discontinuities occur? 

Our study rests on the recent availability of large collections of popular music with asso¬ 
ciated timestamps, and computational methods with which to measure them [25]. Analysis 
in traditional musicology and earlier data-driven ethnomusicology [26], while rich in struc¬ 
ture [18], is slow and prone to inconsistencies and subjectivity. Promising diachronic studies 
on popular music exist, but they either lack scientific rigour [4], focus on technical aspects 
of audio such as loudness, vocabulary statistics and sequential complexity [25,27], or are 
hampered by sample size [28]. The present work uniquely combines the power of a big, 
clearly defined diachronic dataset with the detailed examination of musically meaningful 
audio features. 

To delimit our sample, we focused on songs that appeared in the US Billboard Hot 100 
between 1960 and 2010. We obtained 30-second-long segments of 17,094 songs covering 
86% of the Hot 100, with a small bias towards missing songs in the earlier years. Since 
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Figure 1: Data processing pipeline illustrated with a segment of Queen's Bohemian Rhapsody, 1975, 
one of the few Hot 100 hits to feature an astrophysicist on lead guitar. 
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Figure 2: Evolution of musical Topics in the Billboard Hot 100. Mean Topic frequencies {q) ± 95% 
Cl estimated by bootstrapping. 
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our aim is to investigate the evolution of popular taste, we did not attempt to obtain a 
representative sample of all the songs that were released in the USA in that period of time, 
but just those that were most commercially successful. To analyse the musical properties 
of our songs we adopted an approach inspired by recent advances in text mining (figure 1). 
We began by measuring our songs for a series of quantitative audio features, 12 descriptors 
of tonal content and 14 of timbre (Supplementary Information M.2-3). These were then 
discretised into “words” resulting in a harmonic lexicon (H-lexicon) of chord changes, and a 
timbral lexicon (T-lexicon) of timbre clusters (SI M.4). To relate the T-lexicon to semantic 
labels in plain English we carried out expert annotations (SI M.5). The musical words from 
both lexica were then combined into 8+8=16 “Topics” using Latent Dirichlet Allocation 
(LDA). LDA is a hierarchical generative model of a text-like corpus, in which every document 
(here: song) is represented as a distribution over a number of topics, and every topic is 
represented as a distribution over all possible words (here: chord changes from the H-lexicon, 
and timbre clusters from the T-lexicon). We obtain the most likely model by means of 
probabilistic inference (SI M.6). Each song, then, is represented as a distribution over 8 
harmonic Topics (H-Topics) that capture classes of chord changes (e.g., “dominant 7th chord 
changes”) and 8 timbral Topics (T-Topics) that capture particular timbres (e.g., “drums, 
aggressive, percussive”, “female voice, melodic, vocal”, derived from the expert annotations), 
with Topic proportions q. These Topic frequencies were the basis of our analyses. 


Results 

The Evolution of Topics 

Between 1960 and 2010, the frequencies of the Topics in the Hot 100 varied greatly: some 
Topics became rarer, others became more common, yet others cycled (figure 2). To help 
us interpret these dynamics we made use of associations between the Topics and particular 
artists as well as genre-tags assigned by the listeners of Last.fm, a web-based music discovery 
service with ~50m users (electronic supplementary material, M.8). Considering the H-Topics 
first, the most frequent was H8 (mean ± 95%CI: q = 0.236 ± 0.003)—major chords without 
changes. Nearly two-thirds of our songs show a substantial (> 12.5%) frequency of this Topic, 
particularly those tagged as CLASSIC country, classic ROCK and love (online tables). Its 
presence in the Hot 100 was quite constant, being the most common H-Topic in 43 of 50 
years. 

Other H-Topics were much more dynamic. Between 1960 and 2009 the mean frequency of 
HI declined by about 75%. HI captures the use of dominant-7*'^ chords. Inherently dissonant 
(because of the tritone interval between the third and the minor seventh) these chords are 
commonly used in Jazz to create tensions that are eventually resolved to consonant chords; 
in Blues music, the dissonances are typically not resolved and thus add to the characteristic 
“dirty” colour. Accordingly we find that songs tagged blues or jazz have a high frequency 
of HI; it is especially common in the songs of Blues artists such as B.B. King and Jazz artists 
such as Nat “King” Cole. The decline of this Topic, then, represents the lingering death of 
Jazz and Blues in the Hot 100. 

The remaining H-Topics capture the evolution of other musical styles. H3, for example, 
embraces minor-7th chords used for harmonic colour in funk, disco and soul—this Topic is 
over-represented in funk and DISCO and artists like Chic and KC &: The Sunshine Band. 
Between 1967 and 1977, the mean frequency of H3 more than doubles. H6 combines several 
chord changes that are a mainstay in modal rock tunes and therefore common in artists 
with big-stadium ambitions (e.g., Motley Criie, Van Halen, REO Speedwagon, Queen, Kiss 
and Alice Cooper). Its increase between 1978 and 1985, and subsequent decline in the early 
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1990s, slightly earlier than predicted by the BBC [29], marks the age of Arena Rock. Of ah 
H-Topics, H5 shows the most striking change in frequency. This Topic, which captures the 
absence of an identihable chord structure, barely features in the 1960s and 1970s when, a 
few spoken-word-music collages aside (e.g., those of Dickie Goodman), nearly all songs had 
clearly identihable chords. H5 starts to become more frequent in the late 1980s and then 
rises rapidly to a peak in 1993. This represents the rise of Hip Hop, Rap and related genres, 
as exemplihed by the music of Busta Rhymes, Nas, and Snoop Dog, who all use chords 
particularly rarely (online tables). 

The frequencies of the timbral Topics, too, evolve over time. T3, described as “energetic, 
speech, bright”, shows the same dynamics as H5 and is also associated with the rise of Hip 
Hop-related genres. Several of the other timbral Topics, however, appear to rise and fall 
repeatedly, suggesting recurring fashions in instrumentation. For example, the evolution 
of T4 (“piano, orchestra, harmonic”) appears sinusoidal, suggesting a return in the 2000s 
to timbral qualities prominent in the 1970s. T5 (“guitar, loud, energetic”) underwent two 
full cycles with peaks in 1966 and 1985, heading upward once more in 2009. The second, 
larger, peak coincides with a peak in H6, the chord-changes also associated with stadium rock 
groups such as Motley Crue (online tables). Finally, T1 (“drums, aggressive, percussive”) 
rises continuously until 1990 which coincides with the spread of new percussive technology 
such as drum machines and the gated reverb effect famously used by Phil Collins on In the 
air tonight, 1981. Accordingly, T1 is overrepresented in songs tagged DANCE, DISCO and 
NEW WAVE and artists such as The Pet Shop Boys. After 1990, the frequency of T1 declines: 
the reign of the drum machine was over. 

The varieties of music 

To analyse the evolution of musical variety we began by classifying our songs. Popular music 
is classihed into genres such as country, rock and roll, rhythm and blues (R’n’B) 
as well as a multitude of subgenres (dance-pop, synthpop, heartland rock, roots 
ROCK etc.). Such genres are, however, but imperfect reflections of musical qualities. Popular 
music genres such as country and RAP partially capture musical styles but, besides being 
informal, are also based on non-musical factors such as the age or ethnicity of performers 
(e.g., CLASSIC ROCK and K[orean]-Pop) [5]. For this reason we constructed a taxonomy of 
13 Styles by /c-means clustering on principal components derived from our Topic frequencies 
(hgure 3 and electronic supplementary material M.9). We investigated all k < 25 and found 
that the best clustering solution, as determined by mean silhouette score, was k = 13. 

In order to relate Last.fm tags to the style Style clusters, we used a technique called 
enrichment analysis from bio-informatics. This technique is usually applied to arrive at 
biological interpretations of sets of genes, i.e. to hnd out what the “function” of a set of 
genes is. Applying the GeneMerge enrichment-detection algorithm [30] to our Style data, we 
found that all Styles are strongly enriched for particular tags, i.e. for each Style some Last.fm 
tags are signihcantly over-represented (table SI), so we conclude that they capture at least 
some of the structure of popular music perceived by consumers. The evolutionary dynamics 
of our Styles reflect well-known trends in popular music. For example, the frequency of 
Style 4, strongly enriched for JAZZ, funk, soul and related tags, declines steadily from 1960 
onwards. By contrast. Styles 5 and 13, strongly enriched for ROCK-related tags, fluctuate in 
frequency, while Style 2, enriched for RAP-related tags, is very rare before the mid-1980s but 
then rapidly expands to become the single largest Style for the next thirty years, contracting 
again in the late 2000s. 

What do our Styles represent? Figure 3 shows that Styles and their evolution relate 
to discrete sub-groups of the charts (genres), and hierarchical cluster analysis suggests that 
styles can be grouped into a higher hierarchy. However, we suppose that, unlike organisms of 
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Figure 3: Evolution of musical styles in the Billboard Hot 100. The evolution of 13 Styles, defined by 
fc-means clustering on principal components of Topic frequencies. The width of each spindle is proportional 
to the frequency of that style, normalised to each year. The spindle contours are based on a ±2-year moving 
average smoother; unsmoothed yearly frequencies are shown as grey horizontal lines. A hierarchical cluster 
analysis on the fc-means centroids grouped our Styles into several larger clusters here represented by a tree: 
an EASY-LISTENING -|- LOVE-SONG cluster, a COUNTRY -|- ROCK cluster, and SOUL -I- FUNK -1- DANCE 
cluster; the fourth, most divergent, cluster only contains the hip hop -|- RAP-rich Style 2. All resolved 
nodes have > 75% bootstrap support. Labels list the four most highly over-represented Last.fm user tags 
in each Style according to our enrichment analysis; see electronic supplementary material table SI for full 
results. Shaded regions define eras separated by musical revolutions (figure 5). 


different biological species, all the songs in the charts comprise one large, highly structured, 
meta-population of songs linked by a network of ancestor-descendant relationships arising 
from songwriters imitating their predecessors [31]. Styles and genres, then, represent popu¬ 
lations of music that have evolved unique characters (Topics), or combinations of characters, 
in partial geographic or cultural isolation, e.g., country in the Southern USA during the 
1920s or RAP in the South Bronx of the 1970s. These Styles rise and fall in frequency over 
time in response to the changing tastes of songwriters, musicians and producers, who are in 
turn influenced by the audience. 

Musical diversity has not declined 

Just as paleontologists have discussed the tempo and mode of evolutionary change in the 
fossil record [32], historians of music have discussed musical change and the processes that 
drive it. Some have argued that oligopoly in the media industries has caused a relentless 
decline in cultural diversity of new music [1,2], while others suggest that such homogenizing 
trends are periodically interrupted by small competitors offering novel and varied content 
resulting in “cycles of symbol production” [7,11]. For want of data there have been few tests 
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Figure 4: Evolution of musical diversity in the Billboard Hot 100. We estimate four measures of 
diversity. From left to right: Song number in the charts, Z?n, depends only on the rate of turnover of 
unique entities (songs), and takes no account of their phenotypic similarity. Class diversity, Ds, is the 
effective number of Styles and captures functional diversity. Topic diversity, Dt, is the effective number 
of musical Topics used each year, averaged across the Harmonic and Timbral Topics. Disparity, Dy, or 
phenotypic range is estimated as the total standard deviation within a year. Note that although in ecology 
Ds and Dy are often applied to sets of distinct species or lineages they need not be; our use of them 
implies nothing about the ontological status of our Styles and Topics. For full definitions of the diversity 
measures see electronic supplementary material, M.ll. Shaded regions define eras separated by musical 
revolutions (figure 5). 


of either theory [8-10,13]. 

To test these ideas we estimated four yearly measures of diversity (figure 4). We found 
that although all four evolve, two—Topic diversity and disparity—show the most striking 
changes, both declining to a minimum around 1984, but then rebounding and increasing 
to a maximum in the early 2000s. Since neither of these measures track song number, 
their dynamics cannot be due to varying numbers of songs in the Hot 100; nor, since our 
sampling over 50 years is nearly complete, can they be due to the over-representation of 
recent songs—the so-called “pull of the recent” [33]. Instead, their dynamics are due to 
changes in the frequencies of musical styles. 

The decline in Topic diversity and disparity in the early 1980s is due to a decline of 
timbral rather than harmonic diversity (electronic supplementary material, figure SI). This 
can be seen in the evolution of particular topics (figure 2). In the early 1980s timbral Top¬ 
ics T1 (drums, aggressive, percussive) and T5 (guitar, loud, energetic) become increasingly 
dominant; the subsequent recovery of diversity is due to the relative decrease in frequency 
of the these topics as T3 (energetic, speech, bright) increases. Put in terms of Styles, the 
decline of diversity is due to the dominance of genres such as new wave, disco, hardrock; 
its recovery is due to their waning with the rise of RAP and related genres (figure 2). Con¬ 
trary to current theories of musical evolution, then, we find no evidence for the progressive 
homogenisation of music in the charts and little sign of diversity cycles within the 50 year 
time frame of our study. Instead, the evolution of chart diversity is dominated by historically 
unique events: the rise and fall of particular ways of making music. 

Musical evolution is punctuated by revolutions 

The history of popular music is often seen as a succession of distinct eras, e.g., the “Rock 
Era”, separated by revolutions [3,6,13]. Against this, some scholars have argued that musical 
eras and revolutions are illusory [5]. Even among those who see discontinuities, there is little 
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Figure 5: Musical revolutions in the Billboard Hot 100. A. Quarterly pairwise distance matrix of all 
the songs in the Hot 100. B. rate of stylistic change based on Foote Novelty over successive quarters for 
all windows 1-10 years, inclusive. The rate of musical change—slow-to-fast—is represented by the colour 
gradient blue, green, yellow, red, brown: 1964, 1983, and 1991 are periods of particularly rapid musical 
change. Using a Foote Novelty kernel with a half-width of 3 years results in significant change in these 
periods, with Novelty peaks in 1963-Q4 (P < 0.01), 1982-Q4 (P < 0.01) and 1991-Ql (P < 0.001) 
marked by dashed lines. Significance cutoffs for all windows were empirically determined by random 
permutation of the distance matrix. Significance contour lines with P values are shown in black. 
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agreement about when they occurred. The problem, again, is that data have been scarce 
and objective criteria for deciding what constitutes a break in a historical sequence, scarcer 
yet. 

To determine directly whether rate discontinuities exist we divided the period 1960-2010 
into 200 quarters and used the principal components of the Topic frequencies to estimate a 
pairwise distance matrix between them (figure 5A). This matrix suggested that, while musical 
evolution was ceaseless, there were periods of relative stasis punctuated by periods of rapid 
change. To test this impression we applied a method from Music Information Retrieval, 
Foote Novelty, which estimates the magnitude of change in a distance matrix over a given 
temporal window [34]. The method relies on a kernel matrix with a checkerboard pattern. 
Since a distance matrix exposes just such a checkerboard pattern at change points [34], 
convolving it with the checkerboard kernel along its diagonal directly yields the novelty 
function (SI M.ll). We calculated Foote Novelty for all windows between 1 and 10 years and, 
for each window, determined empirical significance cutoffs based on random permutation of 
the distance matrix. We identified three revolutions: a major one around 1991 and two 
smaller ones around 1964 and 1983 (hgure 5B). From peak to succeeding trough, the rate of 
musical change during these revolutions varied 4- to 6-fold. 

This temporal analysis, when combined with our Style clusters (figure 3), shows how 
musical revolutions are associated with the expansion and contraction of particular musical 
styles. Using quadratic regression models, we identified the Styles that showed significant 
(P < 0.01) change in frequency against time in the six years surrounding each revolu¬ 
tion (electronic supplementary material, table S2). We also carried out a Style-enrichment 
analysis for the same periods (electronic supplementary material, table S2). Of the three 
revolutions 1964 was the most complex, involving the expansion of several Styles—1, 5, 8, 
12 and 13—enriched at the time for SOUL and ROCK-related tags. These gains were bought 
at the expense of Styles 3 and 6 both enriched for doowop among other tags. The 1983 
revolution is associated with an expansion of three Styles—8,11 and 13—here enriched for 
NEW WAVE, DISCO and HARD ROCK-related tags and the contraction of three Styles—3, 7 
and 12—here enriched for soft rock, COUNTRY-related or SOUL -|- R’N’B-related tags. The 
largest revolution of the three, 1991, is associated with the expansion of Style 2, enriched 
for RAP-related tags, at the expense of Styles 5 and 13, here enriched for ROCK-related tags. 
The rise of RAP and related genres appears, then, to be the single most important event that 
has shaped the musical structure of the American charts in the period that we studied. 

The British did not start the American revolution of '64 

Our analysis does not reveal the origins of musical styles; rather, it shows when changes in 
style frequency affect the musical structure of the charts. Bearing this in mind we investigated 
the roles of particular artists in one revolution. On 26 December, 1963, The Beatles released I 
want to hold your hand in the USA. They were swiftly followed by dozens of British acts who, 
over the next few years, flooded the American charts. It is often claimed that this “British 
Invasion” (BI) was responsible for musical changes of the time [35] . Was it? As noted above, 
around 1964 many Styles were changing in frequency; many principal components of the 
Topic frequencies show linear changes in this period too. Inspection of the first four PCs 
shows that their evolutionary trajectories were all established before 1964, implying that, 
while the British may have contributed to this revolution, they could not have been entirely 
responsible for it (figure 6A). We then compared two of the most successful BI acts, The 
Beatles and The Rolling Stones, to the rest of the Hot 100 (figure 6B). In the case of PCI 
and PC2, the songs of both bands have (low) values that anticipate the Hot lOO’s trajectory: 
for these musical attributes they were literally ahead of the curve. In the case of PC3 and 
PC4 their songs resemble the rest of the Hot 100: for these musical attributes they were 
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Figure 6: The British Invasion in the American revolution of 1964. Top to Bottom: PC1-PC4. A. 
Linear evolution of quarterly medians of four PCs in the six years (24 quarters) flanking 1963-Q4, the peak 
of the 1964 revolution. The population medians of all four PCs decrease, and these decreases begin well 
before the start of the British Invasion in late 1963, implying that Bl acts cannot be solely responsible for 
the changes in musical style evident at the time. For each PC, the two topics that load most strongly are 
indicated, with sign of correlation—high, red to low, blue—indicated (electronic supplementary material, 
figure S2). B. Frequency density distributions of four PCs for Beatles, Rolling Stones and songs by all 
other artists around the 1964 revolution. For PCI and PC2, but not PC3 and PC4, The Beatles and 
The Rolling Stones have significantly lower median values than the rest of the population, indicated with 
arrows, implying that these Bl artists adopted a musical style that exaggerated existing trends in the Hot 
100 towards increased use of major chords and decreased use of “bright" speech (PCI) and increased 
guitar-driven aggression and decreased use of mellow vocals (PC2). Vertical lines represent medians; P 
values based on Mann-Whitney-Wilcoxon rank sum test; The Beatles (B): N = 46; The Rolling Stones 
(RS): N = 20; Other artists (0): N = 3114. 
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merely on-trend. Together, these results suggest that, even if the British did not initiate the 
American revolution of 1964, they did exploit it and, to the degree that they were imitated 
by other artists, fanned its flames. Indeed, the extraordinary success of these two groups—66 
Hot 100 hits between them prior to 1968—may be attributable to their having done so. 


Discussion and Conclusions 

Our findings provide a quantitative picture of the evolution of popular music in the USA 
over the course of fifty years. As such, they form the basis for the scientific study of musical 
change. Those who wish to make claims about how and when popular music changed can 
no longer appeal to anecdote, connoisseurship and theory unadorned by data. Similarly, 
recent work has shown that it is possible to identify discrete stylistic changes in the history of 
Western classical music by clustering on motifs extracted from a corpus of written scores [36] . 
Insofar that our approach is based on audio, it can also be applied to music for which no scores 
exist, including that from pre-Modern cultures [18,26,37]. We have already applied a similar 
approach to the classihcation of Art music (“classical music”) into historical periods [38]. 
More generally, music is a natural starting point for the study of stylistic evolution because 
it is not only a universal human cultural trait [39], but also measurable, largely determined 
by form, and available in a relatively standardised format (digital recordings). 

Our study is limited in several ways. First, it is limited by the features studied. Our 
measures must capture only a fraction of the phenotypic complexity of even the simplest 
song; other measures may give different results. However, the hnding that our classifications 
are supported by listener genre-tags gives us some conhdence that we have captured an 
important part of the perceptible variance of our sample. Second, in confining our study to 
the Hot 100, 1960-2010, we have only sampled a small fraction of the new singles released 
in the USA; a complete picture would require compiling a database of several million songs, 
which in itself is a challenge [40]. Given that the Hot 100 is certainly a biased subset of these 
songs, our conclusions cannot be extended to the population of all releases. Finally, we are 
interested in extending the temporal range of our sample to at least the 1940s—if only to 
see whether 1955 was, as many have claimed, the birth date of Rock’n’Roll [41]. 

We have not addressed the causes of the dynamics that we detect. Like any cultural 
artefact—and any living organism—music is the result of a variational-selection process [15- 
18]. In evolutionary biology, causal explanations of organismal diversity appeal to intrinsic 
constraints (developmental or genetic), ecological factors (competition among individuals 
or lineages) and stochastic events (e.g., rocks from space) [42-44]. By analogy, a causal 
account of the evolution of music must ultimately contain an account of how musicians 
imitate, and modify, existing music when creating new songs, that is, an account of the 
mode of inheritance, the production of musical novelty, and its constraints. The first of 
these—inheritance and its constraints—is obscure [45,46]; the second—selection—less so. 
The selective forces acting upon new songs are at least partly captured by their rise and 
fall through the ranks of the charts. Many anecdotal histories of music attempt to explain 
these dynamics. For example, the rise of rap in the charts has been credited to the television 
show Yo, MTV Raps! first broadcast in 1988 [47]. A general, multilevel, selection theory, 
not restricted to Mendelian inheritance, should provide a means for such hypotheses to be 
tested [48-50]. 

Finally, we note that the statistical tools used in this study are quite general. Latent 
Dirichlet Allocation can be used to study the evolving structure of many kinds of assemblages; 
Foote Novelty can be used to detect rate discontinuities in temporal sequences of distances 
based on many kinds of phenotypes. Such tools, and the existence of large digital corpora of 
cultural artefacts—texts, music, images, computer-aided design (CAD) files—now permits 
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the evolutionary analysis of many dimensions of modern culture. We anticipate that the 
study of cultural trends based upon such datasets will soon constrain and inspire theories 
about the evolution of culture just as the fossil record has for the evolution of life [51]. 
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M Materials and Methods 


M.l The origin of the songs 

Metadata on the complete Billboard Hot 100 charts were obtained through the (now defunct) Billboard 
API, consisting of artist name, track name, and chart position in every week of the charts from 1957 to 
early 2010. We use only songs from 1960 through 2009 since these years have complete coverage. Using 
a proprietary matching procedure, we associated Last.fm MP3 audio recordings with the chart entries. 
Each recording is 30 seconds long. We use 17,094 songs, covering 86% of the weekly Billboard charts 
(84% before 2000, 95% from 2000 onward (figure Ml). This amounts to 69% of unique audio recordings. 
The total duration of the music data is 143 hours. 



Figure Ml; Coverage of the Billboard Hot 100 Charts by week. 


To validate our impression that data quality was good, a random sub-sample of 9928 songs was vetted 
by hundreds of volunteers recruited on the internet. The participants were presented with two recordings, 
and for each were asked to to answer the question “Does recording [...] have very poor audio quality?”. We 
analysed those 5593 recordings that were judged at least twice. A recording was considered poor quality 
if it was marked as such by a majority vote. Overall, this was the case in only 3.8% of the recordings, 
with a bias towards worse quality recordings in the 1960s (9.1%; 1970 and later: 1.8%). To examine the 
effect of bad songs, we removed them and compared the estimated mean q of each topic (Section M.6) for 
the total population versus the population of ‘good’ songs for each year of the 1960s. In no case did we 
find that they were significantly different. We conclude that recording quality will have a negligible effect 
on our results. 

All songs were decoded to PCM WAV format (44100Hz, 16bit). The songs were then band-pass- 
filtered using the Audio Degradation Toolbox [1] to reduce differences in recording equalisation in the 
bass and high treble frequencies (stop-band frequencies: 67 Hz, 6000 Hz). 

M.2 Measuring Harmony. 

The harmony features consist of 12-dimensional chroma features (also: pitch class profiles) [2]. Chroma 
is widely used in MIR as a robust feature for chord and key detection [3], audio thumbnailing [4], and 
automatic structural segmentation [5]. In every frame chroma represents the activations (i.e. the strength) 

c = (ci,..., C12) 

corresponding to the 12 pitch classes in the chromatic musical scale (i.e. that of the piano): A, Bb, B C, 
..., G, Ab. We use the NNLS Chroma implementation [6] to extract chroma at the same frame rate as 
the timbre features (step size: 1024 samples = 23ms, i.e. 43 per second), but with the default frame size 
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of 16384 samples. The chroma representation (often called chromagram) of the complete 30 s excerpt of 
“Bohemian Rhapsody” is shown in figure 1 (main text). 

M.3 Measuring Timbre. 

The timbre features consist of 12 Mel-frequency cepstral coefficients (MFCCs), one delta-MFCC value, and 
one Zero-crossing Count (ZCC) feature. MFCCs are spectral-domain audio features for the description 
of timbre and are routinely used in speech recognition [7] and Music Information Retrieval (MIR) tasks 
[8]. For every frame, they provide a low-dimensional parametrisation of the overall shape of the signal’s 
Mel-spectrum, i.e. a spectral representation that takes into account human near-logarithmic perception 
of sound in magnitude (log-magnitude) and frequency (Mel scale). We use the first 12 MFCCs (excluding 
the 0**’ component) and additionally one delta-MFCC, calculated as the difference between any two 
consecutive values of the 0**' MFCC component. The MFCCs were extracted using a plugin from the 
Vamp library (seen 27.03.2014) with the default parameters (block size: 2048 samples = 46ms, step 
size: 1024 samples = 23ms). This amounts to « 43 frames per second. The ZCC (also: zero-crossing 
rate, ZCR) is a time-domain audio feature which has been used in speech recognition [9] and has been 
applied successfully to discern drum sounds [10]. It is calculated by simply counting the number of times 
consecutive samples in a frame are of opposite sign. ZCC is high for noisy signals and transient sounds at 
the onset of consonants and percussive events. To extract the ZCC we also used a Vamp plugin, extracting 
features at the same frame rate (43 per second, step size: 1024 sampes = 23ms), but with a block size 
of 1024 samples. MFCCs and zero crossing counts of “Bohemian Rhapsody” are shown in figure 1 (main 
text). 

M.4 Making musical lexica 

Since we aim to apply topic models to our data (see Section M.6), we need to discretise our raw features 
into musical lexica. We have one timbral lexicon (T-Lexicon) and one harmonic lexicon (H-Lexicon). 

Timbre. 

In order to define the T-Lexicon we followed an unsupervised feature learning approach by qnantising 
the feature space into 35 discrete classes as follows. First, we randomly selected 20 frames from each 
of 11350 randomly selected songs (227 from every year), a total of 227,000 frames. The features were 
then standardised, and de-correlated using principal component analysis (PC A). The PC A components 
were once more standardised. We then applied model-based clustering (Gaussian mixture models, GMM) 
to the standardised de-correlated data, using the built-in Matlab function gmdistribution.fit with full 
covariance matrix [11]. The GMM with 35 mixtures (clusers) was chosen as it minimised the Bayes 
Information Criterion. We then transformed all songs according to the same PC A, scaling and cluster 
mapping transformations. In particular, every audio frame was assigned to its most likely cluster according 
to the GMM. Frames with cluster probabilities of < 0.5 were removed. 

Harmony. 

Our H-Lexicon consists of all 192 possible changes between the most frequently used chord types in 
popular music [12]: major (M), minor (m), dominant 7 (7) and minor 7 chords (m7). We use chord 
changes because they offer a key-independent way of describing the temporal dynamics of harmony. As a 
chord is defined by its root pitch class (A,Bb,B,C,... ,Abt>) and its type, our system gives rise to 4 x 12 = 48 
chords. Each of the chords can be represented as a binary chord template with 12 elements corresponding 
to the twelve pitch classes. For example, the four chords with root A are these. 

CTam = (1,0,0,0,1,0,0,1,0,0,0,0) 
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Figure M2: Chord activation, with the most salient chords at any time highlighted in blue. Excerpt of 
“Bohemian Rhapsody” by Queen. 


CTAm = (1,0,0,1,0,0,0,1,0,0,0,0) 

CTa7 = (1,0,0,0,1,0,0,1,0,0,1,0) 

^,4^7 = (1,0,0,1,0,0,0,1,0,0,1,0) 

At every frame we estimate the locally most likely chord by correlating the chroma vectors (Section ??) 
with the binary chord templates (see, e.g. [13]), i.e. a given chroma frame c = (ci,..., C 12 ) is correlated 
to a chord template CT = {CT \,..., CT 12 ) by using Pearson’s rho, 

^ iCTi-CT){ci-c) 

PCT,c = > --—- 

^ O'CTO'c 

^=l 

where " denotes the sample mean and a. denotes the sample standard deviation of the corresponding 
vector. To smooth these correlation scores over time, we apply a median filter of length 43 (1 second). 
An example of the resulting smoothed chord activation matrix is shown in figure M2. We then choose the 
chord with the highest median-smoothed correlation and combine the two chord labels spaced 1 second 
apart into one chord change label, retaining only the relative root positions of the chords and both chord 
types [14], as demonstrated in figure M2. If the chord change was ambiguous (mean correlation of the 
two chords < 0.4), the chord change label was set to an additional 193’''^ label NA. 

In summary, we have obtained two lexica of frame-wise discrete labels, one for timbre (35 classes) and 
one for harmony (193 chord changes). Each allows us to describe a piece of music as a count vector giving 
counts of timbre classes and chord changes, respectively. 

M.5 Semantic lexicon annotation 

Since we can now express our music in terms of lexica of discrete items, we can attach human-readable 
labels to these items. In the case of the 193 chord changes (H-Lexicon), an intrinsic musical interpretation 
exists. The most frequent chord changes are given in Table ??, along with some explanations and counts 
over the whole corpus. 

The 35 classes in the T-Lexicon do not have a priori interpretations, so we obtained human annotations 
on a subset of our data. First, we randomly selected 100 songs, two from each year, and concatenated 
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the audio that belonged to the same of the 35 sound classes from our T-Lexicon using an overlap-add 
approach. That is, each audio file contained frames from only one of the timbre classes introduced in 
Section M.4, but from up to 100 songs. The resulting 35 sound class files can be accessed on SoundCloud^). 
We noticed that each of the files does indeed have a timbre characteristic; some captured a particular 
vowel sound, others noisy hi-hat and crash cymbal sounds, others again very short, percussive sounds. 
We then asked ten human annotators to individually describe these sounds. Each annotator listened to 
all 35 files and, for each, subjectively chose 5 terms that described the sound from a controlled vocabulary 
consisting of the following 34 labels manually compiled from initial free-vocabulary annotations; 

mellow, aggressive, dark, bright, ealm, energetic, smooth, percussive, guiet, loud, harmonic, noisy, 
melodic, rounded, harsh, vocal, instrumental, speech, instrument: drums, instrument: guitar, instrument: 
piano, instrument: orchestra, instrument: male voice, instrument: female voice, instrument: synthesiser, 
‘ah’, ‘ay’, ‘ee’, ‘er’, ‘oh’, ‘ooh’, ‘sh’, ‘ss’, [random - I find it hard to judge]. 

On average, the most agreed-upon label per class was chosen by 7.5 (mean) of the 10 annotators, 
indicating good agreement. Even the second- and third-ranking labels were chosen by more than half of 
the annotators (means 6.4 and 5.68). Figure M3 shows the agreement of the top labels from rank 1 to 8. 



1 2 3 4 5 6 7 8 

rank of most agreed-upon sound labels 


Figure M3: Agreement of the 10 annotators in the semantic sound annotation task. 

M.6 Topic extraction 

For timbre and harmony separately, a topic model is estimated from the song-wise counts, using the 
implementation of Latent Dirichlet Allocation (LDA) [15] provided in the topicmodels library [16] for 
R. LDA is a hierarchical generative model of a corpus. The original model was formulated in the context 
of a text corpus in which 

a) every document (here: song) is represented as a discrete distribution over Nt topics 

b) every every topic is represented as a discrete distribution over all possible words (here: H-Lexicon 
or T-Lexicon entries) 

Since the T- and H-Lexicon count vectors introduced in Section M.4 are of the same format as word 
counts, we can apply the same modelling procedure. That is, by means of probabilistic inference on the 
model, the LDA method estimates the topic distributions of each song (probabilities of a song using a 
particular topic) and the topics’ lexical distribution (probabilities over the H- and T-lexica) from the 
lexicon count vectors. 

We used the LDA function, which implements the variational expectation-maximization (VEM) algo¬ 
rithm to estimate the parameters, setting the number of topics to 8. Hence, we obtained one model with 

^https;//soundcloud.com/descent-of-pop/sets/cluster-sounds 
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8 T-Topics, and one with 8 H-Topics. Topic models allow us to encode every song as a distribution over 
T- and H-Topics, 


The probabilities can be interpreted as the proportion of frames in the song belonging to the topic. 
When it is clear from the context which T- or H-Topic we are concerned with we denote these by the 
letter q, and their mean over a group of songs by q. Mean values by year for all topics are shown in 
figure ?? in the main text with 95% confidence intervals based on quantile bootstrapping. 

In the same manner, we calculate means and bootstrap confidence intervals for all artists with at least 
10 chart entries and all Last.fm tags (introduced in Section M.8) with at least 200 occurrences. The 
artists with the highest and lowest mean q of each topic and the respective listing of tags can be found 
online. 

M.7 Semantic topic annotations 

In this section we show how to map the semantic interpretations of our harmony and timbre lexica (see 
Section M.5) to the 8 T-Topics and 8 H-Topics. This allows us to work with the topics rather than the 
large number of chord changes and sound classes. 

Harmony. 

Each H-Topic is dehned as a distribution over all H-lexicon entries i = 1,..., 193 (the 193 

different chord changes). The lower half of Table ?? shows the 10 most probable chord changes for 
each topic with those that have P(Ei) > 0.01 emphasised in bold. For example, the most likely chord 
change in H-Topic 4 is a Major chord followed by another Major chord 7 semitones higher, e.g. C to 
G. The interpretation of a topic, then, is the coincidence of such chord changes in a piece of music. 
Interpretations of the 8 H-Topics can be found in Table 1. 

Timbre. 

In order to obtain interpretations for the T-Topics we map the semantic annotations of the T-lexicon 
(Section M.5) to the topics. The semantic annotations of the T-lexicon come as a matrix of counts 
W* = {w*j) of annotation labels j = 1,...,34 for each of the sound classes i = 1,...,35. We first 
normalise the columns w*j by root-mean-square normalisation to obtain a scaled matrix Wij with the 
elements 


Wij — , —^^- 

^(1/34) EiK)" 

The matrix W = {wij) expresses the relevance of the label for the sound class. Since T-Topics are 
compositions of sound classes, we can now simply map these relevance values to the topics by multipli¬ 
cation. The weight Lj of the j**’ label for a T-Topic in which sound class Ef appears with probability 
P{Ef) is 

35 

L,=J2^ijP{E'[). ( 2 ) 

i=l 

The top 3 labels for each T-Topic can be found in Table 1. 
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Table 1: Topic interpretations. 


harmonic topics 

HI changes involving dominant 7th chords 
H2 natural minor key changes 
H3 changes involving minor 7th chords 
H4 simple diatonic changes used in major keys 
H5 unrecognised changes or no chordal content 
H6 stepwise changes indicating modal harmony 
H7 ambiguous major/minor attribution 
H8 sustained major chords 


timbral topics 

T1 

drums, aggressive, percussive 

T2 

calm, quiet, mellow 

T3 

energetic, speech, bright 

T4 

piano, orchestra, harmonic 

T5 

guitar, loud, energetic 

T6 

ay, male voice, vocal 

T7 

oh, rounded, mellow 

T8 

female voice, melodic, vocal 


M.8 User-generated tags 

The Last.fm recordings are also associated with tags, generated by Last.fm users, which we obtained via a 
proprietary process. The tags are usually genre-related (POP, SOUL), but a few also contain information 
about the instrumentation, feel (PIANO, SUMMER), references to particular artists and others. We 
removed references to particular artists and joined some tags that were semantically identical. After the 
procedure we had tags for 16085 (94%) of the songs, with a mean tag count of 2.7 per song (median: 3, 
mode: 4). 

M.9 Identifying musical Styles clusters: A:-means and silhouette scores 

In order to identify musical styles from our data measurements, we first used the 17094 x 16 (i.e. songs 
X topics) data matrix of all topic probabilities and q^, and de-correlated the data using PC A (see 
also figure S2). The resulting data matrix has 14 non-degenerate principal components, which we used 
to cluster the data using fe-means clustering. We chose a cluster number of 13 based on analysing of 
the mean silhouette width [17] over a range of k = 2,... ,25 clusters, each started with 50 random 
initialisations. The result of the best clustering at fc = 13 is chosen, and each song is thus classified to a 
style s 6 {1,..., 13} (figure M4). 
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Figure M4: Mean silhouette scores. The optimal number of clusters, fe = 13 is highlighted in blue. 


M.IO Diversity metrics 

In order assess the diversity of a set of songs (usually the songs having entered the charts in a certain 
year) we calculate four different metrics: number of songs (Htv), effective number of styles (style diversity, 
Dc), effective number of topics (topic diversity, Dt) and disparity (total standard deviation, Dy)- The 
following paragraphs explain these metrics. 

Number of songs. 

The simplest measure of complexity is the number of songs Djy. We use it to show that other diversity 
metrics are not intrinsically linked to this measure. 


Effective number of Styles. 

In the ecology literature, diversity refers to the effective number of species in an ecosystem. Maximum 
diversity is achieved when the species’ frequencies are all equal, i.e. when they are uniformly distributed. 
Likewise, minimum diversity is assumed when all organisms belong to the same species. According to 
[18], diversity for a population of Ng species can formally be dehned as 


Ds = exp 



(3) 


where Sj, i = 1,..., W represents the relative frequency distribution over the Ng species such that Si = 
1. In particular, the maximum value assumed when all species’ relative frequencies are equal is D = Ng. 
If, on the other hand, only one species remains, and all others have frequencies of zero, then D = 1, the 
minimum value. 

We use this exact definition to describe the year-wise diversity of acoustical Style clusters in our 
data (recall that each song has only one Style, but a mixture of Topics). For every year we calculate 
























the proportion songs Si, i = 1,..., 13 belonging to each of the 13 Styles, and hence we use = 13 to 
calculate Ds 6 [1,13]. 


Effective number of Topics. 

The probability g of a certain topic in a song (see Section M.6) provides an estimate of the proportion 
of frames in a song that belong to that topic. By averaging over the year, we can get an estimate of the 
proportion q of frames in the whole year, i.e. for all T- and H-Topics we obtain the yearly measurements 


Figuratively, we throw all audio frames of all songs into one big bucket pertaining to a year, and estimate 
the proportion of each topic in the bucket. From these yearly estimates of topic frequencies we can now 
calculate the effective number of T- and FI-Topics in the same way we calculated the effective number of 
Styles (figure SI). 


D'^ 

Dj' 


= exp 


i=l 
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(4) 

(5) 

( 6 ) 


where we define Dt as the overall measure of topic diversity. Dt is shown in the main manuscript (figure 
4). The individual H- and T-Topic diversities and D!ji are provided in figure ??. ft is evident that 
the significant diversity decline in the 1980s is mainly due to a decline in timbral topic diversity, while 
harmonic diversity shows no sign of sustained decline. 


Disparity. 

In contrast to diversity, disparity corresponds to morphological variety, variety of measurements. Two 
ecosystems of equal diversity can have different disparity, depending on the extent to which the phenotypes 
of species differ. A variety of measures, such as average pairwise character dissimilarity and the total 
variance (sum of univariate variance) [19, 20] have been used to measure disparity. We adopt the square 
root of total variance, a metric called total standard deviation [21, p. 37] as our measure of disparity, i.e. 
given a set of N observations on T traits as a matrix X = "we define it as 


Dy = 


Var(x._„ 




(7) 


We apply our disparity measure Dy to the 14-dimensional matrix of principal components (derived 
from the topics, as described in Section M.9). 


M.ll Identifying musical revolutions 

In order to identify points at which the composition of the charts significantly changes, we employ Foote 
novelty detection [22], a technique often used in MIR [23]. First we pool the 14-dimensional principal 
component data (see Section M.9) into quarters by their first entry into the charts (January-March, April- 
June, July-September, October-December) using the quarterly mean of each principal component. We 
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then construct a matrix (see figure 5 in main text) of pairwise distances between each quarter. Foote’s 
method consists of convolving such a distance matrix with a so-called checkerboard kernel along the main 
diagonal of the matrix. Checkerboard kernels represent the stylised case of homogeneity within regions 
(low values in the upper right and bottom-left quadrants) and dissimilarity between regions (high values 
in the other two quadrants). In such a situation, i.e. when one homogenous era transitions to another, 
the convolution results in high values. 

A kernel with a half-width of 12 quarters (3 years) compares the 3 years prior to the current quarter to 
those following the current quarter (figure M5). We follow Foote [22] in using checkerboard kernels with 
Gaussian taper (standard deviation; 0.4 times the half-width). The kernel matrix entries corresponding 
to the central, “current” quarter are set to zero. 

Many different kernel widths are possible. Figure 5B in the main text shows the novelty score for 
kernels with half-widths between 4 quarters (1 year) and 50 quarters (12.5 years). We can clearly make 
out three major ‘revolutions’ (early 1960s, early 1980s, early 1990s) that result in high novelty scores for 
a wide range of kernel sizes. 



-10 -5 0 5 10 

lag/months 


Figure M5: Foote checkerboard kernel for novelty detection. 

In order to be able to assess the significance of these regions we compared their novelty scores against 
novelty values obtained from randomly permuted distance matrices. We first repeated the process 1000 
times on distance matrices with randomly permuted quarters. For every kernel size we then calculated 
the quantiles at confidence levels a = 0.95, 0.99 and 0.999. The results are shown as contour lines in 
figure 5B in the main text. 

For further analysis we choose the time scale depicted with a half-width lag of 12 quarters (3 years). 
This results in three change regions at confidence p < 0.01 given in Table ??. The ‘revolution’ points are 
the points of maximum Foote novelty within the three regions of significant change, see Table ??. Note 
that there are no significant changes at small time scales (< 2 years). On the other hand, all quarters 
have significant change at large time scales, i.e. the charts evolve long-term. 

M.12 Identifying Styles that change aronnd each revolntion 

To identify the styles (clusters) that change around each revolution, we obtained the frequencies of each 
style for the 24 quarters flanking the peak of a revolution, and estimated the rate of change per annum 
by a quadratic model. We then used a tag-enrichment analysis to identify those tags associated with each 
style just around each revolution, see Table S2. 
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Figure SI: Evolution of Topic diversity. Year-wise topic diversity measures; A D^; B 
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0.0070 


2: calm, quiet, mellow 



a49n 

drums, aggressive, percussive 

1 °D 

1: drums, aggressive, percussive 

-0.22 1 



8: M.O.M, M.5.M, M.9.m7 

-0.43 1 


8: M.O.M, M.5.M, M.9.m7 

-0.021 [ 



7: M.O.m, m.O.M, M.7.m 


0.17 1 

7: M.O.m, m.O.M, M.7.m 

-0.24 I 



6: M.2.M, M.5.M, M.10.M 

-0.21 1 


6: M.2.M, M.5.M, M.10.M 

-«•'» I_ 



5: NA, M.11.M, m.ll.M 


0.43 1 

5: NA, M.ll.M, m.ll.M 

-0.1 1^ 



4; M.7.M, M.5.M, m7.10.M 

-0.32 1 


4; M.7.M, M.5.M, m7.10.M 

-0.091 1 



3: m7.0.m7, M.9.m7, m7.3.M 


(.12 1 

3: m7.0.m7, M.9.m7, m7.3.M 


0.2 1 


2: m.O.m, m.0.m7, M.9.m 


0.17 1 

2: m.O.m, m.0.m7, M.9.m 


0.24 1 


1: 7.0.7, M.0.7, 7.0.M 

0.05 

□ 

1; 7.0.7, M.0.7, 7.0.M 

0.013 




-0.5 0.0 0.5 


-0.5 0.0 

0.5 


Figure S2: Loadings of the first 2 principal components extracted from the Topic probabilities (see 
Section M.9). 
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Table SI: Enrichment analysis: Last.fm user tag over-representation for all Styles over the complete data 
set (only those with P < 0.05) 


Style 1 

Style 2 

Style 3 

Style 4 

Style 5 

Style 6 

Style 7 

northernsoul 

hiphop 

easylistening 

funk 

rock 

femalevocal 

country 

soul 

rap 

country 

blues 

classicrock 

pop 

classiccountry 

hiphop 

gangstarap 

lovesong 

jazz 

pop 

rnb comb 

folk 

dance 

oldschool 

piano 

soul 

newwave 

motown 

rockabilly 

rap 

dirtysouth 

ballad 

instrumental 

garage 

soul 

southernrock 

vocaltrance 

dance 

classiccountry 

rocknroll comb 

hardrock 

glee 


house 

dancehall 

westcoast 

party 

comedy 

reggae 

newjackswing 

classic 

sexy 

urban 

jazz 

doowop 

swing 

lounge 

malevocal 

50s 

singersongwriter 

romantic 

softrock 

Christmas 

acoustic 

northernsoul 

rockabilly 

garagerock 

british 

soundtrack 

musical 


Style 8 

Style 9 

Style 10 

Style 11 

Style 12 

Style 13 


dance 

classicrock 

lovesong 

funk 

soul 

rock 


newwave 

country 

slowjams 

blues 

rnb comb 

hardrock 


pop 

rock 

soul 

dance 

funk 

alternative 


electronic 

singersongwriter 

folk 

bluesrock 

disco 

classicrock 


synthpop 

folkrock 

rnb comb 

newwave 

slowjams 

alt indi rock comb 


freestyle 

folk 

neosoul 

electronic 

dance 

hairmetal 


rock 

pop 

femalevocal 

synthpop 

neosoul 

poppunk 


eurodance 

newjackswing 

triphop 

funk 

soft rock 

singersongwriter 

acoustic 

romantic 

mellow 

easylistening 

jazz 

beautiful 

smooth 

softrock 

ballad 

hiphop 

hardrock 

newjackswing 

smooth 

oldschoolsoul 

punkrock 

punk 

poprock 

metal 

powerpop 

emo 

numetal 

heavymetal 

glamrock 

grunge 

british 
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Table S2: Identifying Styles that change around each revolution and the associated over-represented genre 
tags in the 24 quarters flanking the revolutions. 


revo¬ 

lution 

style 

cluster 

estim. 

(linear) 

P 

estim. 

(quad.) 

P 

tags 

1964 

1 

0.044 

0.019 

0.050 

0.009 

northernsoul, motown, soul, easy listening, jazz 


2 

-0.014 

0.291 

0.008 

0.528 

comedy, funny, jazz, easylistening 


3 

-0.177 

0.000 

0.014 

0.631 

easy listening, jazz, swing, lounge, doowop 


4 

-0.060 

0.118 

-0.009 

0.806 

jazz, northernsoul, soul, instrumental, blues 


5 

0.060 

0.014 

-0.041 

0.081 

garagerock, garage, british, psychedelic 


6 

-0.171 

0.000 

-0.045 

0.282 

femalevocal, motown, northernsoul, soul, doowop 


7 

-0.018 

0.495 

-0.003 

0.908 

british, folk, surf, malevocal, rocknroll comb 


8 

0.062 

0.000 

0.010 

0.441 

garagerock, instrumental, northernsoul, surf, soul 


9 

0.083 

0.014 

0.006 

0.861 

rocknroll comb, northernsoul, motown, soul, garage 


10 

-0.020 

0.349 

0.034 

0.109 

folk, easylistening, jazz, swing 


11 

0.037 

0.131 

-0.022 

0.349 

garagerock, blues, soul, northernsoul, instrumental 


12 

0.043 

0.001 

0.000 

0.996 

northernsoul, soul, motown 


13 

0.121 

0.000 

-0.005 

0.821 

psychedelic, garagerock, psychedelicrock, motown, british 


1983 1 

0.031 

0.350 

-0.050 

0.136 

newwave, synthpop, disco 

2 

0.024 

0.189 

0.010 

0.569 

oldschool, funk, comedy 

3 

-0.112 

0.000 

0.021 

0.426 

lovesong, softrock, easylistening, romantic 

4 

-0.022 

0.448 

-0.005 

0.846 

newwave, disco, progressiverock 

5 

0.020 

0.614 

-0.019 

0.617 

newwave, rock, classicrock, progressiverock, pop, synthpop 

6 

-0.013 

0.598 

-0.005 

0.828 

femalevocal, disco, pop, reggae, musical, soundtrack 

7 

-0.098 

0.003 

0.015 

0.598 

newwave, classiccountry, softrock, newromantic, rock, southernrock 

8 

0.173 

0.000 

-0.029 

0.190 

newwave, pop, rock, disco, synthpop 

9 

-0.065 

0.057 

0.065 

0.055 

classicrock, rock, softrock, progressiverock, newwave, pop 

10 

-0.047 

0.110 

-0.004 

0.882 

lovesong, softrock 

11 

0.048 

0.089 

-0.058 

0.043 

newwave, synthpop, rock, classicrock, hardrock 

12 

-0.095 

0.001 

0.011 

0.664 

funk, disco, soul, smoothjazz, dance 

13 

0.150 

0.000 

0.034 

0.236 

rock, classicrock, hardrock, newwave, progressiverock, southernrock, 
powerpop, heavymetal, hairmetal 


1991 1 

0.034 

0.161 

-0.018 

0.457 

house, freestyle, synthpop, newwave, electronic, dance 

2 

0.325 

0.000 

0.043 

0.187 

hiphop, rap, oldschool, newjackswing, gangstarap, eurodance, westcoast, dance 

3 

0.056 

0.011 

0.011 

0.596 

ballad, lovesong 

4 

0.005 

0.822 

-0.034 

0.153 

dance, hairmetal, hardrock, metal, newjackswing 

5 

-0.085 

0.003 

0.003 

0.901 

rock, hardrock, hairmetal, pop, freestyle 

6 

0.000 

0.999 

0.015 

0.344 

femalevocal, pop, slowjams, dance, rnb comb, ballad 

7 

0.023 

0.307 

0.007 

0.742 

rock, hardrock, hairmetal 

8 

-0.023 

0.584 

-0.023 

0.576 

dance, newjackswing, freestyle, pop, electronic, synthpop, newwave, australian 

9 

-0.015 

0.690 

0.046 

0.224 

rock, pop, softrock, ballad, hardrock 

10 

0.042 

0.235 

0.008 

0.813 

slowjams, lovesong, rnb comb 

11 

-0.065 

0.016 

0.030 

0.248 

newwave, dance, synthpop, hairmetal, rock, freestyle 

12 

-0.035 

0.447 

-0.017 

0.712 

newjackswing, rnb comb, dance, slowjams, house, pop, soul 

13 

-0.207 

0.000 

-0.112 

0.023 

hardrock, hairmetal, rock, classicrock, metal, heavymetal, thrashmetal, 
madchester 


13 





References 


[1] Mauch M, Ewert S. The Audio Degradation Toolbox and its Application to Robustness Evaluation. 
In: Proceedings of the 14th International Society of Music Information Retrieval Conference (ISMIR 
2013); 2013. p. 83-88. 

[2] Fujishima T. Real Time Chord Recognition of Musical Sound: a System using Common Lisp Music. 
In: Proceedings of the International Computer Music Conference (ICMC 1999); 1999. p. 464-467. 

[3] Mauch M, Dixon S. Simultaneous Estimation of Chords and Musical Context from Audio. IEEE 
Transactions on Audio, Speech, and Language Processing. 2010;18(6):1280-1289. 

[4] Bartsch MA, Wakefield GH. Audio thumbnailing of popular music using chroma-based representa¬ 
tions. IEEE Transactions on Multimedia. 2005;7(1):96-104. 

[5] Muller M, Kurth F. Towards structural analysis of audio recordings in the presence of musical 
variations. EURASIP Journal on Applied Signal Processing. 2007;2007(1):163-163. 

[6] Mauch M, Dixon S. Approximate note transcription for the improved identification of difficult chords. 
Proceedings of the 11th International Society for Music Information Retrieval Conference (ISMIR 
2010). 2010;p. 135-140. 

[7] Davis S, Mermelstein P. Comparison of parametric representations for monosyllabic word recognition 
in continuously spoken sentences. IEEE Transactions on Acoustics, Speech and Signal Processing. 
1980;28(4):357-366. 

[8] Foote JT. Content-Based Retrieval of Music and Audio. Proceedings of SPIE. 1997;138:138-147. 

[9] Ito M, Donaldson RW. Zero-crossing measurements for analysis and recognition of speech sounds. 
IEEE Transactions on Audio and Electroacoustics. 1971;19(3):235-242. 

[10] Gouyon F, Pachet F, Delerue O. On the use of zero-crossing rate for an application of classification of 
percussive sounds. In: Proceedings of the COST G-6 Conference on Digital Audio Effects (DAFX-00); 

2000. p. 1-6. 

[11] McLachlan GJ, Peel D. Finite Mixture Models. Hoboken, NJ: John Wiley & Sons; 2000. 

[12] Burgoyne JA, Wild J, Fujinaga 1. An expert ground-truth set for audio chord recognition and 
music analysis. In: Proceedings of the 12th International Conference on Music Information Retrieval 
(ISMIR 2011); 2011. p. 633-638. 

[13] Papadopoulos H, Peeters G. Large-scale Study of Chord Estimation Algorithms Based on Chroma 
Representation and HMM. In: International Workshop on Content-Based Multimedia Indexing; 
2007. p. 53-60. 

[14] Mauch M, Dixon S, Harte C, Fields B, Casey M. Discovering chord idioms through Beatles and Real 
Book songs. In: Proceedings of the 8th International Conference on Music Information Retrieval 
(ISMIR 2007); 2007. p. 255-258. 

[15] Blei DM, Ng AY, Jordan ML Latent dirichlet allocation. Journal of Machine Learning Research. 
2003;3:993-1022. 

[16] Hornik K, Griin B. topicmodels: An R package for fitting topic models. Journal of Statistical 
Software. 2011;40(13):l-30. 

[17] Rousseeuw PJ. Silhouettes: a graphical aid to the interpretation and validation of cluster analysis. 
Journal of computational and applied mathematics. 1987;20:53-65. 


14 



[18] Jost L. Entropy and diversity. Oikos. 2006; 113(2):363-375. 

[19] Erwin DH. DISPARITY: MORPHOLOGICAL PATTERN AND DEVELOPMENTAL CONTEXT. 
Palaeontology. 2007;50(l):57-73. 

[20] Foote M. The evolution of morphological diversity. Annual Review of Ecology and Systematics. 
1997;. 

[21] Hallgrimsson B, Hall BK. Variation: A Central Concept in Biology. Academic Press; 2011. 

[22] Foote J. Automatic audio segmentation using a measure of audio novelty. In: IEEE International 
Conference on Multimedia and Expo. vol. 1; 2000. p. 452-455. 

[23] Smith JBL, Chew E. A Meta-Analysis of the MIREX Structure Segmentation Task. In: Proceedings 
of the 14th International Society for Music Information Retrieval Conference; 2013. . 


15 



